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SUMMARY V 

Calculations have "been, made for the icing limit of a. diamond air- 
foil at zero angle of attack in 1 terms of the stream Mach number, stream 
temperature, and pressure altitude. The icing limi t i-s- defined as a 
wetted-surface temperature of 32° F and is related to the stream condi- 
tions "by the method of Hardy. . 

The results show that the point most likely to . ice on the airfoil 
lies immediately behind the shQulder and is subject j to possible icing 
at Mach numbers as high. as 1.4« 

INTRODUCTION . ' 

As the operational speed of aircraft, is increased through the 
transonic region, the frictional heating available to prevent the forma- 
tion of ice on the aircraft becomes an important quantity. A study of 
the probable icing limits of a' high-speed airfoil in transonic and super- 
sonic flow was made in terms of the pertinent variables of flight Mach 
number and free- stream conditions' of pressure and temperature. The pro- 
cedure used to determine the icing limit is baked oh the method presented 
in reference 1 for calculating the surface, temperature of an insulated 
body (a body within which heat is not conducted from one section to 
another) running fully wet in an air stream. The iCing limit at any 
point on the airfoil was assumed to depend upon maintaining a tempera- 
ture of 32° F at that point. The free-stream static rtemperature corre- 
sponding to this icing limit was calculated for each point on the airfoil 
for particular values of free -stream Mach number and pressure. 

The results presented herein were calculated for a symmetrical 
diamond airfoil at zero angle of attack for a range of airfoil- 
thickness ratios from 0.02 to 0.10, pressure altitude from sea level to 
45,000 feet, and free -stream static temperatures to -40° F. 
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SYMBOIB 

The following symbols are used in this report: 
c-n specific heat of air at constant pressure, Btu/(lh) (°F) 

sr * 

e vapor pressure, Ib/sq ft 
kg coefficient of evaporation 

k^ coefficient of heat transfer 

L latent heat of vaporization, Btu/lb 

M Mach number 

* 

molecular weight of air 
m e molecular weight of water vapor 

p static pressure, Xb/sq ft 

r recovery factor 

T static temperature, °R 

t/c airfoil thickness ratio, ratio of airfoil maximum thickness to 
chord 

Tq c mini mum free -stream static temperature corresponding to ice -free 
condition on surface as defined in equation (la), °R 

AT-^ temperature increment resulting from frictional heating in the 
Boundary layer, °Rr 

X ratio of specific heats of air, 1-400 
Subscripts : 

0 free -stream static conditions 

1 local conditions at the edge of Boundary layer 

s surface 


2704 



2704 


MCA Tiff 2861 


3 


AMLYSIS 

The method of reference 1 offers a convenient me ans for determining 
the local surface temperature on a wetted’ body as a function of the 
free-stream conditions of temperature, pressure, vapor pressure, and 
velocity, provided the flow field about the body is known. Ease of cal- 
culation was achieved by a number of simplifyin g assumptions- The 
restrictions caused by these assumptions will be discussed in relation 
to the general heat balance for an unheated body flying through icing 
conditions . 

Several investigators have made quite complete analyses of the heat 
balance for a body flying through a cloud (references 2 and 3) • The 
generalized heat balance for such conditions is well known and has been 
stated quite concisely by Messinger (reference 2) . Some simplification 
of the general equation may be obtained by considering the case of an 
insulated body subjected to sufficient frictional heating to maintain 
the surface ice free. In this case, the heat balance is as follows: 

(l) Heat due 'to the frictional, or viscous, effects 

plus 

(2) Heat due to the kinetic energy of the water droplets 

plus 

(3) Heat of fusion 
equals 

(4) Heat lost by convection 
plus 

(5) Heat for evaporation of water 
plus 

(6) Heat required to raise temperature of water droplets ' from stream 
temperature to surface temperature 

The result of Hardy's relation (reference l) is obtained by equat- 
ing the frictional term (l) to the sum of the convective term (4) and 
the evaporation term (5). For Hardy's relation to be valid, terms (2) 
and (3) must be small compared with term (l), and term (6) must be small 
compared with term (4) or (5); or the sum of terms (2) and (3) must be 
nearly equal to term ( 6) • 
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The heht of fusion term (3) offers considerable complication to the 
problem. It is conceivable that the' .surface film might consist of a 
mixture of -water and ice particles without ice actually adhering to the 
surface- The heat, qf, fusion would necessarily have to be withdrawn from 
a portion of -the film 'for such pdrticles to occur and, because of the 
lack of knowledge of the percentage of the .film that might be frozen 
before visible icing, occurs, it would appear advisable to define an ice- 
free surface as being -fully wet 'at 32° F with no ice particles in the 
filni. This definition for the icing' limit eliminates term (3) from the 
general heat balances and ‘Hardy 's relation becomes applicable if term (2) 
is approximately ’equal to -term (6) -or both are negligible. Term (2) is 
directly proportional to the impingement- rate and the square of the 
flight speed; term ■(&)• is directly proportional to the impingement rate . 
and the temperature, difference' between ‘the stream and the surface. If 
the impingement rate is small, 'as for 'the case of low liquid -water 
contents, then botja ‘terms (2)- and (6) . are negligible. For calculations 
of the icing limit, both terms . are of the 'same order of magnitude and 
because of the nature. 6f the problem," the temperature difference between 
the surface and. the. stream increases with increasing flight speed- Both 
terms therefore .increase with increasing Mach number; large values for 
either term are only -eneo.Unt eted 'at large values of liquid-water content 
and flight speedy Such a 'combination is highly improbable in actual 
flight since high speeds usually 'occur at high altitudes, and high alti- 
tudes imply small liquid- water contents. 

The relation of 'JEfardy applies in particular to a surface barely 
wetted with a yery thin - film free of ice particles but is probably 
nearly correct for the whole range of practical interest for high-speed 
flight. . 


Method of Calculation 

Wetted-surface temperature . - The equation for the temperature of 
a wetted surface in ah air stream developed in reference 1 was written 
as , , - . 


. T-l - T p + AE-l = 


k e m a / . e s 
k hFaW " fs 



( 1 ) 


where the subscript s refers to conditions at the surface and sub- 
script 1 to conditions at the edge of the boundary layer. From the 
energy relation for a viscous fluid 


T 1 + hT 1 = 



r (r - 1) 
2 



Substituting this relation in equation (l) and solving for T s yields 
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T = 


- L 1 + 


Hr. - i) 


% 




-J hPafy \Pi - e s Pi ' e l' 


Because interest in this problem is' fixed on a definition of flight 
circumstances that provide local surface temperatures of 32° F, the 
terms T s , L, and e s are constants.. Hardy observed (reference l) that 
for the range of temperatures hear . 32° F the ratio of the evaporation 
coefficient k e to the heat -transfer coefficient kh was very nearly 
1. A solution of equation (l) can be obtained for any one of the four 
variables (eq, pq, Mq, and Tq) provided the 'other three are known. 

From purely physical reasoning, however, an additional restriction exists. 
This restriction results from the fact that the air in a cloud is fully 
saturated at the static or free-Btream conditions. For each value of 
stream-static temperature, the fr.ee -stream vapor pressure is therefore 
assumed a constant and equal to the' saturated vapor pressure. If it is 
assumed that the flow about the body, outside, the boundary layer, is 
accomplished with no change in phase, that is, no condensation or evap- 
oration, then Dalton's law of partial pressures applies and 

fl 1 ' '. 

■ e o i>o ■ ■ 


If a reasonable value of the recovery factor r is assumed, say 
0.88, then equation (l) maybe rewritten in the form shown below 


492 = 


Tq (.1 + 0.176 Mq 2 ) 


% 


0,c 


T 0 , c - 2776,5 


12.75 


i 0 


Pq t 12 • 75 


(Po " e 0^ 


(la) 


which is quite convenient for purposes, of calculation because 
T^/Tq^ c (l + 0.176 Mq 2 ) and ViJvq car "be readily calculated for each 

point on a diamond airfoil. 

Because eQ is a function of. To, c> the solution of equation (la) 
must be solved by trial arid error:. In an effort to reduce the laborious - 
ness of such calculations, a number of charts similar to that shown in 
figure 1 were made. For each pressure altitude considered, the free- 
stream static temperature correspon din g to an ice-free condition T q^ c 

was plotted as a function of the parameter Tq/To,c ( 1 + 0.176 Mq 2 ) by 
means of equation (la) for constant values of the pressure ratio Pq/Po' 
For each value of pq/p 0 and Tq/TQ^ C (l + 0.176 Mq 2 ) a free-stream 

static temperature which corresponds fro a surface temperature of 32° F 
and therefore represents the minimum free-stream temperature for an ice- 
free surface Tq, c is obtained. The lower limit of free-stream static 
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temperature was considered to "be -40° F for the purposes of these cal- 
culations since it has "been shown (references 4 and 5) that supercooled 
water droplets are not lik ely to he found at temperatures "below this 
value . 

Flow-field calculations . - The flow field about a diamond or double- 
wedge profile can be characterized in general by four regions; a subsonic 
region, a subsonic region with sonic velocity at the shoulder, a mixed 
subsonic-supersonic region with a detached bow wave, and a supersonic 
region with attached shock waves. The regions of most interest are in 
the currently important ranges of transonic and supersonic flight and 
the results presented herein will be limited thereto. 

A recent analysis is given in references 6 and 7 for the sub- 
sonic region with a sonic line at the shoulder. This solution is for 
a thin, symmetrical wedge at zero angle of attack followed by a straight 
section. The results in this region presented herein are therefore 
limited to the front half of the diamond airfoil. 

The mixed subsonic- supersonic region with detached bow wave has been 
solved in references 8 and 9. This solution applies to a thin, symmet- 
rical diamond airfoil at zero angle of attack for Mach numbers from 1 to 
the Mach number at which the shock attaches. The flow in this region is 
subsonic over the front half of the airfoil increasing to sonic velocity 
at the shoulder. The flow then expands supersonically around the corner, 
and supersonic flow exists over the back half of the airfoil to the 
trailing edge. 

The flow about a diamond airfoil in purely supersonic flow has been 
known for many years- The results used herein are those of reference 10 
because of their ease in use and applicability. The flow in this region 
is supersonic over the front half of the airfoil with an expansion 
region at the shoulder and a consequent higher Mach number on the back 
surface of the airfoil, the velocities and pressures being constant over 
each surface for a particular value of airfoil thickness and stream Mach 
number. 

The calculation of the chordwise pressure and Mach number distri- 
butions over a symmetrical-diamond airfoil of specific thickness ratio 
t/c at zero angle of attack through the entire transonic- supersonic range 
can therefore be accomplished by combining the results of references 6 
to 10. 

The analyses of references 7 to 9 were made using a speed' function £ 
and a generalized pressure coefficient Cp. The results of these analyses 
are presented graphically by plots of the chordwise distribution of £ 
and Cp for constant values of the free -stream speed function £q. The 
pressure ratio P^/Pq is related to the generalized pressure coefficient 

Cp by the equation 
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Pi 

— = i + 

*>o 

2 (r + i) 3 

The Mach number, either local or free stream, is related to the 
speed function £ by the equation 


M = \ 1 + £ 


<r.+ «(D] 


The temperature ratio T^/Tq is obtained from the energy relations 

for a compressible fluid and is given in terms of Mach number by the 
equation 


T 1 + 
1 • 


T, 


0 1 


r - i 

2 

r - i 


Mo 


M, 


An example- of a typical calculation' of Tq^ c at a particular 

chordwise position is given for the 25-percent chard station and a value 
of the free-stream speed function £ 0 of 0.703, Values of £ of -0-740 
and Cp of 2.890 are obtained from figures 3 and 5, respectively, of 
reference 8. If t/c of 0-06 is assumed, then Mq = 1-092, Mi = 0.8926, 
P-j/Pq = 1.276, Tp/To » 1.068, and' "Tj/Tq (l + 0.176 Mi 2 ) = 1-218. At a 

pressure altitude of 15,000 feet, a T 0 c of 422.5° R is obtained from 
figure 1- ' * 

For the supersonic case -with an attached shock (reference 10) , the 
values of Pp/Po & ad M^/Mq can be obtained directly from the charts of 

the reference report and, hence, the calculation of the parameters 
required for obtaining Tq^ c from figure 1 is simple. 


RESULTS, AIJ® DISCUSSION 


For each value of airfoil thickness ratio and free-stream Mach 
number assumed, the values' -of the parameters T-j/Tq (l + 0-176 Mq 2 ) and 

Pl/p 0 can be calculated over the -whole range of transonic and supersonic 


speeds. At each pressure altitude^ "the fre^-> stream static temperature 
Tq^ c corresponding to a wetted surface temperature of 32° F (Ice free) 

can then be obtained from charts similar to figure. 1- 
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This value of Tq^ c corresponds to the lowest free-stream static 

temperature Tq which will give ice -free conditions at the particular 
point on the airfoil under consideration for the given values of free- 
stream Mach number and pressure altitude. 

A plot of T 0 ^ c - as a function of the chordwise distance is shown 

in figure 2 for several free-stream Mach numbers for a 4-percent-thick 
symmetrical-diamond airfoil at zero angle of attack and a pressure alti- 
tude of 30,000 feet. 

A study of figure 2 shows that the. general trend of the curves is 
very similar for Mach numbers from 0. 7593 : to 1.105, that is, from the 
subsonic region to the supersonic region with attached shock. For each 
curve the minimum Tq c increases slowly from- the leading edge to the 

midchord. The sudden jump in the free-stream static temperature required 
for icing protection at the midchord results from the expansion around 
the corner and the consequent decrease in local static pressure and 
increase in' local velocity. The increase in velocity and the decrease 
in pressure enhance the evaporation, hence, the increase in the required 
To, c - Water impinging on the front half of the airfoil might therefore 
run hack along the surface, spill around the shoulder, and freeze. 
Downstream of the shoulder, the free-stream temperature required for 
protection decreases slowly to the trailing- edge. 

When, the shock wave is attached and, the flow is everywhere super- 
sonic (Mq = 1.2, or greater), the free-stream temperature for protection 
is a constant for the whole front surface and the whole hack surface. 

The free-stream temperature required for protection of the hack half is, 
however, greater for the reasons previously cited. 

This and similar figures for other altitudes and airfoil thickness 
ratios show that the critical points for these airfoils are immediately 
ahead of the shoulder for the front half of the airfoil and immediately 
behind the shoulder for the hack half of the airfoil. The icing l imi t 
for these two critical points is shown in figure 3(a) where the minimum 
free-stream temperature corresponding to an ice -free condition is shown 
as a function of the free-stream Mach number for pressure altitudes of 
sea level and 30,000 feet for a 6 -percent -thick airfoil. The value of 
Tq^ c decreases rapidly with increasing Mach number,' that is, the tem- 
perature region of possible icing hazard becomes much smaller with 
increasing Mach number. The free-stream static temperatures required 
for protection apparently are higher at the higher pressure altitudes as 
would he expected since, with other conditions equal, the decreased 
pressure enhances the evaporation rate. At a pressure altitude of 
30,000 feet the curve shown for the critical point behind the shoulder 
indicates the possibility of icing at a Mach number approaching 1.4. 
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The effect of altitude is shown clearly in .figure 3(h) for free -stream 
Mach numbers of 0 - 67 , 1.0, and 1.3 for the critical point ahead of the 
shoulder on the front surface and for Mach numbers of 1.0 and 1.3 for 
the point imme diately behind the shoulder. At a particular value of 
Mach number the curves for both critical points remain nearly parallel 
over the altitude range. For Mach numbers of 1.0 or less the free- 
stream temperature required for protection increases with increasing 
altitude with a definite decrease in the slope of the '-curves in the 
higher altitude range- For the higher Mach numbers the slope of the 
curves increase over the whole range from sea' level to 45,000 feet. 

The effect of airfoil thickness on the icing limit for the two 
critical points on the airfoil is shown in figure. 3(c) where the free- 
stream temperature required for protection is plotted as a function of 
the airfoil thickness ratio. .' The icing limit at a pressure altitude of 
45,000 feet is shown 'for the point ahead of the shoulder for free-stream 
Mach numbers. of 0-875, 1.00, and 1-30 and for the point immediately 
behind the shoulder for Mach numbers of 1.0 and 1.3- For the point 
ahead of the shoulder, the temperature required for protection remains 
constant or decreases with increasing thickness depending on the Mach 
number. For the point immediately behind the shoulder, the temperature 
required for protection increases with increasing thickness. 

An increase in Tq c of approximately 20° R resulted when the 

thickness ratio increased from 0-02 to 0*10 for the critical point 
behind the shoulder at a free-stream Mach -number , of 1-00. • A decrease in 
Tq c of approximately 14.5° R resulted when the thickness ratio 

increased from 0.02 to 0.10 for the critical point ahead of the shoulder 
at a free-stream Mach number of 1.30. 

Experimental verification of the results presented herein has been 
obtained from data taken in a 4- by- 10 -inch duct tunnel at a free-stream 
Mach number of 1.36 .with a 6. 82-percent-ihick symmetrical-diamond airfoil 
of 5.50-inch chord at zero angle of attack. The icing limit was deter- 
mined by successively raising the free-stream static temperature until 
no ice formed on the surface. The results of these tests are summarized 
in the following table and compared with results calculated by the 
analytical method presented herein: 


Free-stream 

static 

pressure 

*o 

(in. Hg) 

Minimum free-stream ; 
corresponding to ice 

static temperature 
-free condition, Tq^ c 

Experimental 

Analytical 

Front 

surface 

Back 

surface 

Front 

surface 

Back 

surface 

9.13 

404 



403 

— 

9.41 

— 

422 

— 

422 

9.76 

406 

— 

401 - 

— 

9.91 



418 

— 

422 

3.91 



— 

447 

— 
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The excellent agreement obtained in several cases is probably somewhat 
fort ui tous since the experimental error was estimated to be ±10 F. 

These results, however, show that the calculated icing limit is very 
closely approximated' experimentally and, hence, provides considerable 
credence for the calculations presented herein. The run-bach icing on 
the rear surface, which was I n d i cated by the analysis , was also observed 
experimentally. Visual observations of the model showed that the surface 
water film flowed smoothly around the shoulder and quickly froze. The 
results presented in the table show that the free-stream static tempera- 
ture at which’ the run -back freezing occurred exper l m entally is in good 
agreement ■with the analytical values. 

The potential hazard to aircraft and missile operation in the tran- 
sonic region may be seen by applying the results previously given for 
the critical points at the shoulder of a diamond airfoil to operation 
under atmospheric conditions likely to be encountered" If a particular 
relation of pressure and temperature with altitude in the atmosphere is 
assumed, it is possible to calculate the flight Mach number necessary 
for icing protection for' a given airfoil thickness as a function of 
altitude. The Mach number necessary for protection for the critical 
point ahead of the shoulder for a 6 -percent-thick diamond airfoil is 
shown in figure 4 for MCA standard atmosphere (reference ll) . For these 
conditions the region of icing hazard exists from approximately 7500 to 
27,300 feet. If a more severe temperature relation with altitude is used 
such as might he encountered during winter, that is, a 32° F sea -level 
surface tenperature ■ and a lapse rate of 3° F per thousand feet, then the 
curves shown result for the critical points at the shoulder. For these 
conditions the region of hazard exists from sea level to 24,000 feet. 

A Mq of 1.29 is needed at the peak altitude in order to provide 

protection. 

Altho ugh the results presented in this report apply only to the 
symmetrical-diamond airfoil at zero angle of attack, they may he used 
as an indication of the region of icing hazard and hence serve as a 
guide for determining where de-icing systems are needed or what flight 
plan may he used to prevent icing. 


SUMMARY OF RESULTS 

The results of the analysis of icing limit presented for the 
diamond airfoil at zero angle of attack in the transonic -supersonic 
speed range are as follows: 

1. The critical point for the whole airfoil in terms of the icing 
limit was located immediately behind the shoulder. The critical point 
on the front surface was located immediately ahead of the shoulder. 

The results indicate that icing may occur at the critical point behind 
the shoulder for Mach numbers as high as 1.4. 
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2. At a constant free -stream Mach number and airfoil thickness, the 
minimum free-stream static temperature corresponding to an ice-free 
condition increases -with increasing altitude, 

3. For the point immediately ahead of the shoulder the minimum 
free-stream temperature corresponding to an ice-free condition remains 
the same or decreases with increasing thickness ratio at a constant 
value of Mach number and altitude. For the point immediately down- 
stream of the shoulder the free-stream temperature corresponding to 

an ice-free condition increases with increasing thickness at a constant 
Mach number and altitude. 


Lewis Flight Propulsion Laboratory 

National Advisory C ommi ttee for Aeronautics 
Cleveland, Ohio, July 24, 1952 
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Figure 2. - ChordHise variation of free-stream fltatic temperature corresponding to ioe-free condition 
for several Mach numbers. Airfoil thiokness, 4 peroentj preoBure altitude, 30,000 feet. 
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Pressure altitude, ft 

(b) Effect of pressure altitude for several Mach numbers, Airfoil thickness, 6 percent. 

Figure 3. - Continued. Variation of minimum free-stream static temperature corresponding to ice- free condition 

for tvo critical noin+s on air foil - 
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Figure 4. - Variation of flight iJach number required for icing protection as 
function of altitude for several assumed atmospheric temperature-pressure 
relations. Airfoil thickness,, 6 percent. 
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